Negative scaling dimensions and conformal invariance at the Nishimori point in the 

±J random-bond Ising model 
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We reexamine the disorder-dominated multicritical point of the two-dimensional ± J-Ising model, 
known as the Nishimori point (NP). At the NP we investigate numerically and analytically the 
behavior of the disorder correlator, familiar from the self-dual description of the pure critical point 
of the two-dimensional Ising model. We consider the logarithmic average and the qth moments of 
this correlator in the ensemble average over randomness, for continuous q in the range < q < 2.5, 
and demonstrate their conformal invariance. At the NP we find, in contrast to the self-dual pure 
critical point, that the disorder correlators exhibit multi-scaling in q which is different from that of 
spin-spin correlators and that their scaling dimension becomes negative for q > 1 and q < 0. Using 
properties on the Nishimori line we show that the first moment (q = 1) of the disorder correlator is 
exactly one for all separations. The spectrum of scaling dimensions at the NP is not parabolic in q. 



PACS numbers: 05.70.Jk, 75.10.Nr, 75.40.Mg 
I. INTRODUCTION 

Conformal field theory (CFT) has been a key in- 
gredient in understanding critical models and phase- 
transitions in two spatial dimensions^. It has helped 
the exact solution of a number of problems and made 
it possible to construct a classification of critical points. 
While so far the success of CFT has been mainly in 
its application to pure critical systems, one might hope 
that it will be an equally powerful tool for understand- 
ing disorder-dominated critical points in two dimensions, 
such as those in disordered magnets or the integer quan- 
tum Hall effect: It is generally expected that the observ- 
ables of a random system at criticality exhibit conformal 
invariance when averaged over disorder, and the averages 
are believed to be vacuum expectation values of opera- 
tors in an appropriate conformal field theory. Even if this 
theory is unknown, the bare assumption of the existence 
of an underlying conformal field theory leads to strong 
constraints which may be tested for particular models. 

One model system which has received renewed atten- 
tion recently is the ±J random-bond Ising model (±J 
RBIM), which (amongst others) has a special symmetry 
line in parameter space, known as the Nishimori linen. It 
has a disorder-dominated multicritical point, the Nishi- 
mori point (NP), where this line intersects the phase 
boundary between a fcrromagnct and a paramagnet. The 
NP may be a good candidate for the construction of a 
consistent CFT of a random critical point since it is one 
of the simplest model critical points at the outset. The 
phase transition in the RBIM is also of interest as a ver- 
sion of the quantum Hall plateau transition, since there 
is an exact mappingdH from the RBIM to a jaetwork 
model similar that used to represent the latteJa. Fur- 
thermore, owing to extensive numerical and analytic ef- 
fort invested in studying the ±J RBIM and the NP in 
particular, accurate estimates for various exponents are 
available, which can be used to test candidate CFTs. At 
the NP the first few integer q moments of the spin-spin 



(crcr-) correlation function have been calculated numeri- 
cally recentlyu'H and have been demonstrated to obey the 
conformal constraints very accurately. The g-dependence 
of their scaling dimension has considerable significance, 
and one can imagine three possibilities. If critical be- 
havior at the NP were that of the percolation transition, 
the dimensions would be independent of q. By contrast, 
the critical behavior in a pure system gives dimensions 
increasing linearly in q. In fact, at the NP scaling dimen- 
sions of moments of the era- correlator vary non-linearly 
with q. This behavior is generic in the presence of ran- 
domness and referred to as multi-scaling. For the RBIM, 
Read and LudwigB considered the expectation values of 
the disorder operator which is the Kramers- Wannier dual 
to the conventional order operator a. The dual operator 
will be the main focus of this paper and we will refer 
to it as ^-operator throughout in order to avoid confu- 
sion between thermal and configurational disorder. The 
/i-opcrators are associated with plaquettcs of the Ising 
lattice: the two-point correlator (^(x)fi(y)) in a particu- 
lar distribution of bonds can be represented by choosing 
an arbitrary path connecting the pair of plaquettes at 
x and y and reversing the sign of all the bonds crossed 
by the path. In this way one arrives at a modified sys- 
tem with partition function Z' (compared to the Z of 
the unmodified mocLcl) and the the correlator is given by 
the ratio of the two!, (fx(x)fx(y)) = Z'/Z. In Ref. | the 
authors point out that the /^-operator, is, in contrast to 
a, not bounded from above anymore when antiferromag- 
netic exchange interaction are present. This opens up 
the interesting possibility of correlation functions which 
increase with distance and consequently the possibility of 
negative scaling dimensions, which in turn is intimately 
linked to the central problem of non-unitarity in CFT. 
The fact that /i/i-corrclators can increase with distance 
has been established in Ref. || for RBIM's with equal con- 
centrations of ferromagnetic and antiferromagnetic ex- 
change interactions. Here our concern is instead with 
behaviour at a critical point. The aim of this paper is 
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to explicitly demonstrate conformal invariance with neg- 
ative scaling dimensions by considering the /u/z-correlator 
at the NP. One practical problem in these calculations is 
to achieve high sampling while keeping the system size 
sufficiently large. From this perspective making use of 
the recently developed mappingu of the RBIM onto the 
network model for calculating free energies is crucial since 
it reduces a formerly exponentially large calculation to 
one which is only power-law in system size. 

In addition to our numerical results, we show in this 
paper that the distribution functions of the /i-operators 
satisfy strong constraints on the Nishimori line. In par- 
ticular we consider the qth moments of correlation func- 
tions and we show that the moments are symmetric in q 
with respect to q = 1/2. Owing to this symmetry the [i- 
operators for q = 1 average exactly to unity, as discussed 
in the following section. Due to the asymmetry between 
the spin operator and the disorder operator in the pres- 
ence of bond disorder, the NP is not microscopically self- 
dual. Despite this the question arises whether self-duality 
is restored in the correlation functions asymptotically for 
large separations. The above exact result on the Nishi- 
mori line shows that this does not happen, since the spin- 
spin correlation function decays with distance. In Sec. Ill 
we present results from extensive numerical calculations 
on the correlations of the disorder /i-operator for the ± J 
RBIM. We show that the moments of its correlation func- 
tion obey the conformal constraints thus reinforcing the 
idea of an underlying conformal field theory. We find 
also that the gth moment of the /z/i-correlation function 
increases with separation r for q > 1, establishing the 
presence negative scaling dimensions and non-unitarity. 
We test our results against previous calculations in the 
quasi one-dimensional regime and find excellent agree- 
ment. Finally, our results show that multi-scaling at the 
NP is not parabolic. 



II. DISTRIBUTIONS IN THE RBIM ON THE 
NISHIMORI LINE 

We consider the two-dimensional nearest-neighbor 
Ising model on the square lattice with partition function 
Z({ J}, 0) = Tr CT exp[J2^j) 0Jij<?i(Jj], where a t is a classi- 
cal spin variable taking the values ±1 and the (ij) denote 
nearest neighbors. The exchange couplings Jy are drawn 
independently from a probability distribution P(J). For 
general P(J) this model is known as the random-bond 
Ising model (RBIM). In the conventional notation the 
correlator of a local operator O x is 

(O x O y ) = |Tr CT O x O y exp[J2 0Jij^jl (2.1) 



It was shown by Kadanoff and Cevacl, however, that if O 
represents either the order operator or its dual operator 
an alternative way of writing Eq.(|2.lD is to absorb the 



product O x O v into the Hamiltonian by modifying the 
set {J^} — * {Jij} so that 



(O x O y ) = Z'/Z, 



(2.2) 



where Z' is the partition function evaluated with the set 
{Jij}- One important observation when considering av- 
erages over bond distributions is that the properties of 
the average of (O x O y ) crucially depend on whether the 
modification O x O y induces takes the set {Jij} out of the 
random ensemble or not. To be specific, inserting the 
order operator a of the Ising model may be viewed as 
giving the bonds along a semi-infinite path on the lat- 
tice an imaginary part while the /i-operator is equivalent 
to simply reversing the bond signs along a semi-infinite 
pathEl. For the /i-order operator the modified set {J[j} is 
still within the physical ensemble. Its average is there- 
fore bounded below by zero but it can take any positive 
valuoJ, in contrast the order parameter which lies be- 
tween — 1 and +1. 

Let us now consider the average of an observable 
W({ J}, 0). W might represent a two-point correlator as 
in Eq.(2T) but is generally an arbitrary function over the 
bond configurations {J} and of temperature. Ensemble 
averages will be denoted by an over-bar 



W = 



W({J},0)HP(Jij)dJi. 



(2.3) 



It was noticed by NishimoriB that this average may be 
rewritten as 

W = Tr T / W({rJ}, 0)Z({0~ 1 tA}, 0) \{ Q(J t] )dJ l3 
Jo ij 

(2.4) 

where Z({(3~ 1 tA}, 0) is an Ising partition function with 
a set of couplings {0~ 1 TijA i3 } and 

Q{Jij) = [P(Jij) + P{-Jij)]/{2co8hAijl (2.5) 
A(Jij) = ~ ]n[P(Jij)/P(-Jij)}, (2.6) 

where ry = ±1. Thus the integral is taken over the pos- 
itive values of Jy only and the negative part of P(Jy) 
is included through the sum over ry. The Nishimori 
line is defined as the line on which A(Jij) = 0Jij, and 
hence Z^^tA}, 0) = Z({tJ},0). Owing to this re- 
lation a number of averages can be performed exactly 
for bimodal and Gaussian distributions (and, others) on 
the Nishimori line. In particular NishimoriB calculated 
the internal energy in these two distributions exactly by 
considering the average of the free energy, In Z. The sym- 
metry on the Nishimori line may also be expressed as a 
replica symmetrytJ from which further strong constraints 
on the distribution functions of correlators follow. 

Now, for the qth moment of a disorder-depen dent ob- 
servable such as the correlation function in Eq.(2.2) on 
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the Nishimori line one can write 

(O x O y )i = TY T ^ i^—j ZJlQiJijWij. (2.7) 

From this it is evident that there appear symmetries in 
the distribution of such correlation functions if there is 
a simple relation between Z and Z' . In particular since 
any product of /i-operators is represented by reversing 
the sign of a set of bonds, Eq.(2.7) is invariant under 
the change Z — > Z 1 ' . This can be seen by noticing that 
taking the trace over r inside the integral takes positive 
and negative bond strengths with equal weight. In turn, 
the change Z — > Z' is equivalent to q — > 1 — q and hence 
the moments are symmetric with respect to q = 1/2; 
The result (OiOj) q — 1 for q = 0, 1 readily follows from 
normalization of Eq.(2.4). 



where 



III. NISHIMORI POINT OF THE ±J-RBIM 

The result of the previous section holds for the gen- 
eral RBIM with a Nishimori line. Here we restrict our- 
selves to the ± J model with the bimodal bond distribu- 
tion P(Jij) = pS(Jij — 1) + (1 — p)5{Jij + 1) for which 
the Nishimori line is given by exp(— 2/3) = (1 —p)/p. On 
this line, the model has a disorder-dominated multicriti- 
cal point known as the Nishimori point. Its position has 
been estimated in the past and again, more accurately, 
recentlyu'EI. Two critical exponents are now known from 
numerical calculations!! We also focus on the Nishimori 
point here and we consider the gth moments of the /Lt/x- 
correlator with separation r. We denote the averaged 
correlators by g(r; q). If they are critical with power law 
decay in the plane then 



g{r;q) = (wrY = A{r/a)-^\ 



(3.1) 



where A and a are constants and 77(g) is a (/-dependent 
scaling dimension. We would like to answer the following 
questions: Is this averaged correlator conformally invari- 
ant? If so, how are the {fJ,ofJ> r ) distributed, or, equiva- 
lently, what is the functional dependence of 77(g) on g? 
As shown in the previous section g(r; q) — 1 for all r at 
q = and q = 1. The scaling dimension 77(g) of g has 
then two zeros in q. Furthermore 77(g) must be less than 
or equal to zero for g > 1 (and g < 0) as shown by the fol- 
lowing argument: we start from the inequality x p > (x) p 
for any real random x satisfying < x < 00 and any real 
p > 1. Applying this to x = (/UoMr) one readily finds that 
v(p) — PV(1) an d, since 77(1) = 0, the above result fol- 
lows. Below we use numerical calculations to show that 
77(g) is, in fact, non-zero and negative for q > 1 and we 
determined its form in detail for a range of q. 

Conformal fiel d th eory predicts that the correlation 
functions in Eq.(3.1) decay on a cylinder with circum- 
ference M according to 



f(r/M) = sm(irx/M) or sinh(7ry/M) 



(3.3) 



and r = (x,y) with x and y the coordinate separations 
around the circumference and along the cylinder, respec- 
tively. In addition one can also consider the logarithmic 
average (typical value) which is related to the spectrum 
in q t hrou gh the derivative with respect to g at g = 0. 
If Eq.(3__2) holds then the typical value decays on the 
cylinder as 



ln(/J /-O 



( drjjg) 
[ 8q 



ln[(M/o7r)/(r/M)]. (3.4) 



Using the mapping of the RBIM onto a network model 
of free fermions described in detail in Ref. ^, and ex- 
ploiting the efficiency of the transfer-matrix algorithm in 
the network form alism we calcu late the correlation func- 
tions in Eq.( |3.2| ) and Eq.(3.4). The logarithmic aver- 



( M 0Mr) 9 = [(M/an)f(r/M)r^ 



(3.2) 



age has small fluctuations compared to the direct aver- 
ages and can be used as a sensitive test of whether the 
two-point correlators are actually conformally invariant. 
Fig. ([l]) shows the dependence of the typical value of the 
/^-correlator on cylinders with circumferences between 
M = 6 and M = 22 with separation along the cylinder 
(x = 0), calculated using 10 4 samples. The statistical 
errors here are between 0.5 and 1.0 percent for all sys- 
tem sizes and hence the error bars in the plot appear 
smaller than the symbols. The data fall accurately onto 
a single straight line with slope k = —0.6911(17). The 
vertical dashed line is the point r = M where, roughly, 
the correlator crosses over from two-dimensional to quasi 
ID behavior. The typical correlation function hence con- 
vincingly obeys the prediction from conformal invariance. 

The direct averages in Eq.(|3.2|) have extremely large 
fluctuations and are dominated by rare events. There- 
fore, compared to computing the logarithmic average 
many more samples are needed. In the following we have 
typically used 10 8 disorder realizations. We calculate the 
average correlation functions on a cylinder now with the 
separation r running around the circumference (y = 0) 
for q- values between g = and g = 2.5; Four of the di- 
rect averages for q = 0.40,0.60,1.00,1.28 and M = 22 
are displayed in Fig.(|J). Again, error bars are smaller 
than the symbol sizes but now there is a strong system 
size dependence: the statistical errors range from 0.5 per- 
cent for the smallest q to 4.0 percent for q = 1.28. The 
functional form dictated by Eq.(|3.2j) is, again, obeyed 
accurately and 77(g) can be read off as the slope of the 
fitted lines. In Fig.(^|) the slopes for various values of q 
are plotted versus g. As expected the curve is symmetric 
with respect to g = 1/2 and has a zero at g = 1, crossing 
to negative values for q > 1. The fact that we reproduce 
77(1) = accurately provides a strong check of our nu- 
merical calculations, since this is a property that emerges 
only after averaging over random bond realisations. 

Our error margin for 77(g) shown in Fig. (^J) is less than 
three times the symbol sizes for q > 1 but smaller than 
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FIG. 1. Scaling of the typical averages of the /j/j-correlation 
functions for M = 6 — 22 where a straight line represents the 
form predicted by conformal invariance. The vertical dashed 
line marks the crossover point from two-dimensional to quasi 
one- dimensional behavior. The statistical errors are between 
0.5 and 1.0 percent throughout. 



this for q < 1. Since the direct averages are dominated 
by rare events, under-sampling must be a concern. For 
higher q this practical problem becomes worse. In order 
to assess sample-to-sample fluctuations we have therefore 
repeated the data analysis while omitting different sub- 
sets of the data: through this we arrive at the above error 
margin. 

Differentiating the curve 77(g) numerically at q — 
in order to compare the result with the typical average 
gives drj/dq » 0.70(1) which agrees with the above value 
of k. On a cylinder of circumference M this part of the 
multi-scaling spectrum may also be compared to previ- 
ous results in the limit of infinite separation. For the 
/i/x-correlator in the Ising model this limit amounts to 
the insertion of an infinite seam of reversed bonds along 
the cylinder. Then, since the logarithmic average of Z 
is simply the free energy, the typical value of the infi- 
nite r /i/i-correlator Z'/Z is the free energy difference 
between periodic and anti-periodic boundary conditions. 
The free energy per site A/ is the interfacial tension 
which conformal field theory predicts to be airM~ 2 with 
a a universal scaling amplitude. For the ± J RBIM at the 
Nishimori p^dntit was shown by independent numerical 
calculations^^ that Af w 0.691(2)7rM~ 2 . This agrees 
within errors with the derivative at q = and coincides 
very well with k = —0.6911(17) obtained above. 

In the context of multifractality it is interesting to 
know whether the spectrum in Fig.(||) is piecewise ex- 
actly parabolic. Such behavior is known to occur for 
the inverse participation ratios of the ground state wave= 
functions of fermions in a random gauge potentialtj 
and it has been proposed and confirmed numerically for 
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wave-functions at the integer quantum Hall transitionEj. 
Parabolic spectra stem, in two dimensions, from quan- 
tities whose logarithms have the distribution of a free 
massless field, and it can be shown that this would be 
compatible with the constraints on the Nishimori line. 
However, while for q < 1.2 the g-dependence of r\ on q 
seems indeed compatible with a parabola, at larger q the 
curve 77(g) versus q deviates significantly and an exactly 
parabolic spectrum at the Nishimori point can thus be 
excluded. This is shown in Fig.(|^) where the dashed line 
is a parabola fitted to 77(g) between q = and q = 1. 



IV. CONCLUSIONS 

We have demonstrated the presence of negative scaling 
dimensions at the multi-critical point in the ± J random- 
bond Ising model by considering the the Kramers- 
Wannier dual operator. Since the latter it is not bounded 
from above in the presence of antiferromagnetic bonds its 
correlation functions can increase with separation in this 
case which we have demonstrated numerically for the ± J 
random-bond Ising model in two dimensions. Further- 
more at the disorder dominated multi-critical point in 
this model we find that the gth moments of the two-point 
correlator in the ensemble distribution behave according 
to the prediction from conformal invariance and that, for 
g > 1, the exponents 77(g) become negative. We have cal- 
culated the multiscaling spectrum 77(g) and we support 
our observations using symmetry properties of the Nishi- 
mori line on which the critical point is located and which 
forces 77(0) = 77(1) = 0. 
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FIG. 3. The scaling dimension of the g-moment of the 
/u/n-correlator as a function of q. The dashed line shows the 
approximate parabolic spectrum valid for 0.2 < q < 1.2. For 
larger q the deviation is significant. 
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